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Perturbative estimates of operator coefficients for improved lattice actions are becoming increasingly important 
for precision simulations of many hadronic observables. Following previous work by Dimm, Lepage, and Mackenzie, 
we consider the feasibility of computing operator coefficients from numerical simulations deep in the perturbative 
region of lattice theories. Here we introduce a background field technique that may allow for the computation 
of the coefficients of clover-field operators in a variety of theories. This method is tested by calculations of 
the renormalized quark mass in lattice NRQCD, and of the 0(a 3 ) clover coefficient for Sheikholeslami-Wohlert 
fermions. First results for the coefficient of the magnetic moment operator in NRQCD are also presented. 



Recent simulations of many hadronic observ- 
ables using tadpole-improved actions, on both 
coarse and fine lattices, have yielded promising 
results. This has spurred further development of 
ever more highly improved lattice actions. At the 
same time, it has also become apparent that one 
must go beyond tree-level improvement in order 
to obtain precision results for many quantities. 

The increasing diversity and complexity of 
improved actions makes traditional perturba- 
tive calculations of operator coefficients in these 
actions problematic. A simpler alternative is 
needed, if only to provide reasonable estimates of 
these coefficients. Dimm, Lepage, and Mackenzie 
have shown that the renormalized mass for Wil- 
son fermions, and the energy zero in NRQCD, can 
be calculated from numerical simulations, deep in 
the perturbative region of these theories Q. 

In this work we introduce a background field 
technique for perturbative matching of lattice 
theories by numerical simulation. We use this 
technique to estimate the coefficients of clover- 
field operators for both Sheikholeslami-Wohlert 
§ and NRQCD g fermions. 

The clover coefficients are tuned by computing 
quark propagators in a uniform background mag- 
netic field, and matching the spin-flip energy 

AE(gB) = E(l,gB)-E(tgB) (1) 



in the lattice theory with the same quantity in 
continuum QCD. The continuum calculation has 
been done by Sapirstein Q , using the background 
field formalism || . We used a lattice formulation 
of the background field method. This method has 
the important feature that the combination gB 
of bare coupling g and bare field strength B does 
not get renormalized. Thus lattice and continuum 
results for AE can be directly compared, without 
regard to the different regulators that are used. 

We used the Wilson gluon action, and split 
the link variable into a classical background part 
and a fluctuating dynamical part, Utot,/j,{x) = 
Udyn,/j.(x)U c \,[i(x) . For an Abelian background 
(as used in the continuum U c ^ ll —2{x\) = 

e igBx!\ 3 _ p er i oc ii c boundary conditions also re- 
quire Ud^ix! = N,X 2 ) = e -i9BNx 2 X 3j and 

gB = 27rn cxt /A 2 , where n cxt is an integer. 

This background field is an exact solution of the 
equations of motion for the effective action, with 
zero external current. This is because the only co- 
variant vector operators that can be formed from 
the classical field strength F c \„ v are all identically 
zero (e.g., D ^[A^F^^ = 0). Consequently we 
do not need to make a perturbative calculation of 
the external current required to support a more 
general classical field. 

On the other hand, a given background field 
does not minimize the effective action. In a nu- 
merical simulation one observes "tunneling" of 



2 



the system, from the extremum at -F c i, to the 
minimum at zero field. To compute observables 
in the presence of the classical field we modify the 
path integral, introducing a "current" J into the 
action as an intermediate step, in order to create 
a minimum in the effective potential near F c \: 



Z[J} = J[dU dyn 



e -P(S wa [U dyxl U i]+\jJ[U dyn ]) _ ( 2 ) 



Physical quantities in the presence of the classi- 
cal field are obtained by extrapolation in the cou- 
pling, Aj — > 0. However, one must perform this 
extrapolation from sufficiently large Aj, in order 
for the system to be dominated by fluctuations in 
the dynamical fields near the classical field. 

To test this procedure, the NRQCD spin 
flip energy was computed twice, using two cur- 
rents, Ji = J2 X tj. Re - TT i u <iyn,n(x)}, and J 2 = 
(V^)E^Im-Tr[^dyn, M .Fci,i2]- The extrapo- 
lated results agree, as shown in Fig. [l]. 
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Figure f . NRQCD spin-flip energy for two exter- 
nal currents (f3 — 16, Mq = 1.8). 

The action (and hence the current J) must be 
invariant under background field transformations, 

U dy n ifl (x) -> fl(x)U dyn ft\x), ( 3 ) 

in order to preserve the nonrenormalization of 
gB. Likewise, we require a background field- 
covariant gauge fixing prescription for the dynam- 
ical fields, when computing quark propagators: 

{0\T ty(x)$(0)) |0) = 

J [dU^e-^+^K-^Ug^x] ,(4) 



where U dyn is the gauge-fixed dynamical field. We 
gauge fix with respect to transformations of the 
dynamical field (with fixed classical field) 

u dyn ^(x) ^ n(x)u dyn u cl (x)n\x + fip^x). (5) 

We used a background field axial (U dyn ^ = I) 
plus Coulomb C=i D[U c \]iU dyn s = 0) gauge fix- 
ing pL These lattice transformations and gauge 
fixings have the desired continuum limits Q (note 
that the quantum symmetry, Eq. (|5|), is broken 
by the external current just as in the continuum). 

We used this approach to calculate the renor- 
malized quark mass M rcn in lattice NRQCD. Rel- 
ativistic corrections up to 0(v A ) in the mean 
quark velocity || were included. The renormal- 
ized masses were extracted from quark correlation 
functions with small nonzero three-momentum. 
Sample results are shown in Fig. ||. All simu- 
lations were done on 12 3 x 24 lattices at p = 9, 
with n cxt = 2, and current J\ was used. Sim- 
ulations were done both with and without tad- 
pole improvement (using the average plaquette 
to compute uq in the former case). 
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Figure 2. NRQCD renormalized quark mass and 
spin-flip energy vs. Aj (P = 9, M = 1.18, with 
tadpole improvement). 

We did linear extrapolations to Aj = to ob- 
tain the "true" masses, given in Tables | and |. 
Errors include a systematic error due to the ex- 
trapolation. The Monte Carlo ( "MC" ) estimates 
of the masses are in good agreement with results 
from perturbation theory ("PT"), due to Morn- 
ingstar jfj). Notice that most of the renormaliza- 



tion is due to tadpoles, which is clearly resolved 
by the Monte Carlo results. 

Table 1 

Renormalized quark mass and clover coefficient in 
lattice NRQCD for three bare masses Mo, without 
tadpole improvement. 
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Simulation results for AE in NRQCD arc also 
shown in Fig. ^. We find the coefficient of the 
magnetic moment operator, C4, by matching the 
(extrapolated) lattice AE to the continuum value 
Q. Defining C4 = 1 + c^a s , and using a s = 
ay(7r/a), we obtain the results in Tables [l]and ^. 

To calibrate our NRQCD calculations, we 
analyzed spin-flip energies for Sheikholeslami- 
Wohlert (SW) fermions, where the clover coef- 
ficient c sw is known 0. We used a tree-level 
clover coefficient, and the value determined non- 
perturbatively for massless fermions in Ref. 
(c sw — 1.24 at (3 — 9). We worked at modest 
quark masses (met « 0.5) at which 0((ma) 2 ) er- 
rors may dominate over 0(a s ma) errors. This 
prevents a direct matching of the lattice results 
with the continuum. Instead, we compare the 
simulation results for AE with the spin-flip en- 
ergy for "classical" SW fermions [Ud yn = I, 
c S w — 1, and with the classical mass identified 



with E(1 , gB, Aj) measured in the simulation]. 
We expect the effects of 0(a 2 ) errors will roughly 
cancel in this approximate matching procedure. 

Sample results are shown in Fig. ^[ The sim- 
ulation results with a tree-level clover coefficient 
show a significant departure from the matching 
condition. At this quark mass (met w .5), a value 
of c sw somewhat smaller than the coefficient for 
massless fermions is suggested by these results. 
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Figure 3. Monte Carlo AE for SW fermions, 
compared to the "classical" spin-flip energy AE C \ 
{fl = 9, k = .115). 

In summary, we have introduced a background 
field method for simulations in the perturbative 
region of lattice theories. Encouraging prelimi- 
nary results were obtained for the coefficients of 
clover operators for SW and NRQCD fermions, 
and for the renormalized quark mass in NRQCD. 
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